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Abstract — Spectrum sensing is a fundamental problem in 
cognitive radio. We propose a function of covariance matrix 
based detection algorithm for spectrum sensing in cognitive 
radio network. Monotonically increasing property of function of 
matrix involving trace operation is utilized as the cornerstone 
for this algorithm. The advantage of proposed algorithm is 
it works under extremely low signal-to-noise ratio, like lower 
than -30 dB with limited sample data. Theoretical analysis of 
threshold setting for the algorithm is discussed. A performance 
comparison between the proposed algorithm and other state-of- 
the-art methods is provided, by the simulation on captured digital 
television (DTV) signal. 

Index Terms — Function of matrix, spectrum sensing, low SNR 
detection, cognitive radio network 

I. Introduction 

AS a limited nature resource, wireless spectrum becomes 
increasingly scarce due to the evolution of various wire- 
less technique. However, it is not utilized efficiently, the 
current utilization of a licensed spectrum varies from 15% to 
85% H]. And the number is even lower in rural area. Cognitive 
radio (CR) is a key technology to mitigate the overcrowding 
of spectrum space based on its capability of dynamic spectrum 
access (DSA). 

One fundamental requirement of CR system is to find the 
frequency hole by secondary user (SU). Therefore, spectrum 
sensing plays an important role in CR system. Each SU should 
be able to sense primary user's (PU) existence accurately in 
low signal-to-noise ratio (SNR) to avoid interference. In some 
sense, spectrum sensing can be treated as signal detection 
problem. Generally, matched filter detection, energy detection, 
cyclostationary feature detection, and covariance matrix based 
detection jS), O, ID, Q are used for spectrum sensing. 
Covariance matrix related spectrum sensing algorithms were 
extended by employing multiple antennas at the cognitive 
receiver Q. Feature template matching (FTM) I?] extracts 
signal feature as the leading eigenvector of signal's covari- 
ance matrix. Kernel principal component analysis (PCA) and 
kernel generalized likelihood ratio test (GLRT) [8| based on 
matched subspace model map original space data into a higher 
dimensional feature space. 
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However, it is difficult to solve the sensing problem under 
extremely low SNR, like -30 dB. To circumvent this difficulty, 
function of matrix based detection (FMD) algorithm will be 
employed first time for spectrum sensing, which does not 
require the prior knowledge of structure of either original 
signal or noise. If PU signal is included in the received signal, 
the trace operation on function of received covariance matrix 
will always be larger than the function of pure noise covariance 
matrix. Thus we can detect the existence of PU. 

II. System Model 

In this paper, we consider there is one receive antenna to 
detect one PU signal. Let x{t) be the continuous-time received 
signal after unknown channel. Let Tg — 1//^ be the sampling 
period, the received signal sample is x[n] = x {nTg). There 
are two hypotheses to detect PU signal's existence, T-Lq, signal 
does not exist; and "Hi, signal exists. The received signal 
samples under the two hypotheses are given respectively as 
follows: 

T-Lq :x [n] =w [n] (1) 
Hi -.x [n] = s[n]+ w [n] (2) 

where w [n] is the received white noise, and each sample 
of w [n] is assumed to be independent identical distribution 
(i.i.d.), with zero mean and variance a^^. s [n] is the received 
PU signal samples after unknown channel with unknown 
signal distribution. And, it is assumed that noise and signal 
are uncorrected. 

Assume spectrum sensing is performed based on the statis- 
tics of the i*'* sensing segment Tx,i, which consists of Ns 
sensing vectors with L (called "smoothing factor") consecutive 
output samples in each vector: 

^x,i = {x(i_i)Ar^ + i, X(i_i)7v^_(_2, • • • ,^{i-l)N^+N^^ (3) 

= [x[i\ ,x[i + I] , ■ ■ ■ ,x[i + L - l]f (4) 

where ~ M {Q,'R.x)- If Ns is large enough, R^; can be 
approximated by sample covariance matrix R^;: 

^- = ]^E^'^^ (5) 

We will use Rj. instead of R^. for convenience. Accordingly, 
we have R^ for S; and R„ for Wj. As we know, R^ is a 
positive semi-definite matrix with low rank, while R„ is a 
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positive definite matrix. If received signal contains original 
signal and noise, then we have: 

R:r = + R„ + E (sw^) + E (ws"^) (6) 

Since the data samples employed is large enough, we assume 
signal and noise to be uncorrected, which implies E (^sw'^) — 
E (ws'^^ — 0, then 

R-T = R-s + R-n (7) 

So our two hypotheses can be rewritten as 

'Hq : Rx — Rn (8) 

Hi : lUx = Rs + Rra (9) 

The detection will be based upon the covariance matrix of 
received signal in sensing segment T^^i- 

III. Function of Matrix Based Detection 

A. Spectrum Sensing Algorithm 

The total received sensing segment T^^i is divided into K 
sub-segments, with Nk vectors in each sub-segment, where 
Nk — Ns/K. Tx^i can be represented as 

^x,i = Wl,'P2, ■ ■ ■ ,IPk} (10) 

The /c*'* sub-segment 7^ is represented as 

fk = {x(fe_i)Ar^ + i,X(fc_i)jVfc+2, • • • ,^{k~l}Nk+Nk} (H) 

Then we calculate the sample covariance matrix from each 
sensing sub-segment separately, using Eq. (jSj) with sum inter- 
val from 1 to TVs replaced by different sub-segment interval. 
Covariance matrix of fc*'* sub-segment is denoted as 



Rt 



{k^l)Nk+Nk 

Nk 



(12) 



i={k-l)Nk + l 

Accordingly, we have original signal covariance matrix R^.fe 
and noise covariance matrix R„fc. Similar to Eq. Q, if 
original signal is contained, we have 

R£c,fc = Rs,fe + Rn,fc (13) 

Because Hg^k is positive semi-definite matrix, for all k 

Rs,fc + Rn.fc >~ Rn.fc (14) 

In order to make the metric robust, we calculate the sample 
mean of the covariance matrix to get a average covariance 
matrix. From Eq. ([14]), we get 

1 ^ 

, (Rs,fc + Rn,/c) >- JT ^ Rn,fe (15) 
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Generally, if A B, we may find many appropriate 
monotonically increasing functions applied to both sides of the 
matrix inequality with matrix inequality holds that / (A) -< 
/ (B), as function of matrix is still matrix. But not all the 
monotonically increasing functions satisfy / (A) -< f (B). 
Fortunately, if we apply trace operation to both sides of the 
function /, the real number inequality always holds, as shown 
in Fact 1. p. 529]. 



Fact 1. Let A,B e F"""" (F" means R" or C"), assume 
that A and B are positive semidefinite, assume that A ^ B, 
let / : [0, oo) t-^ [0, oo), and assume that /(O) = 0, / is 
continuous, and / is increasing. Then, 



Tr(/(A))<Tr(/(B)) 



(16) 



Based on Eq. (15 i and ([T6]l, we can get 

Tr(^/(^^g(Rs,fc + R„,fc))) 



>Tr(/( igR„.fc 



(17) 



We can see from the above results that if PU signal exists, the 



left side of Eq. (17 i will be larger than right side. Because we 
have averaged the covariance matrices of sub-segments, the 
numerical difference between them is stable. As a result, we 
can find a threshold to distinguish the PU signal present or 
absent. 

Based on the analysis above, we propose detection algo- 
rithm 1 for spectrum sensing. Namely, FMD. 

Algorithm 1 FMD 

1: Determine probability of false alarm Pfa and threshold 7. 

2: Divide the total received segment into K sub-segments, 
and calculate their covariance matrices. 

3: Get an average covariance matrix from all covariance 
matrices of sub-segments 

4: Find an appropriate rational function / and get the func- 
tion / result of average covariance matrix. 

5: Obtain the trace of the function of covariance matrix as a 
metric p 



Tr 



(18) 



if p > 7 then 

PU signal exists 

else 

PU signal does not exist 
end if 



B. Detection and Threshold Setting 

The optimal choice of appropriate monotonically increasing 
function is data dependent. For simplicity, we choose function 
/(A) = A, which is monotonically increasing and can achieve 
near best performance among all the potential functions we 



experienced. Thus, the metric Eq. (I81 can be specified to 



/ 1 

\ k=l 



(19) 



In Eq. ^T9\ , all the covariance matrices are added together 
firstly, then Tr operation is applied. Since Tr operation is noth- 
ing but a sum operation of diagonal elements of covariance 
matrix. We can apply Tr operation to covariance matrix firstly. 
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then add them together to get the same pi, which can be 
represented as 



1 ^ 

Pi = -r^^^'^{^x,k) 



(20) 



k = l 



In this paper, we use Eq. ( 20 1 to determine probability of false 
alarm and the threshold. 

For hypothesis Ho, original signal is not contained in 
received signal. Since the number of vectors to calculate 
covariance matrix is very large, based on Central Limit 
Theorem, ^Tr (R^;) can be approximated by Gaussian distri- 
bution with mean and variance 



151. Then 



or 2 



Pi 



k = l 



N Lai 



2I?_ 



(21) 
(22) 



The variance of noise ct„ can be estimated from Eq. (21 



k 



2L2 



(T^ is the sample variance of Tr (R^, fc), which can be obtained 



from unbiased estimator: 



K 



—.tiZk-ZK) 



k=l 



Where, 



Zk — Tr (Rx^k) 
1 ^ 

Zk = -^Tr(R,,fc) 



fc=i 



(24) 

(25) 
(26) 

The choice of threshold 7 is a compromise between Pd and 
Pfa- Usually, it is difficult to set threshold based on Pd, so we 
choose the threshold based on Pfa- The probability of false 
alarm for FMD is 



Then a„ is replaced by (7„ in the Eq. (22 1 to be 



P. 



fo. 



Pfa = p (pi > 7 1 no) 

p I Pl - L(fn^ > " ^'^"^ 



2L^ - 4 



2L^ 4 



where 



Pfa = 




7 - Ld„ 



2L^ ' 4 



+ 00 



Then we get the threshold 7 

l={l + Q-'{Pfa) 



2 da:: 



KNk 



(28) 
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IV. Covariance Matrix Based Algorithms 

Several sample covariance matrix based algorithms have 
been proposed in spectrum sensing. Maximum-minimum 
eigenvalue (MME) Q and arithmetic -to-geometric mean 
(AGM) |4| have no prior knowledge, while FTM Q has 
feature as prior knowledge. Estimator-Correlator (EC) requires 
both original signal covariance matrix and noise variance. All 
the threshold 7 are determined by probability of false alarm 
if not specified. 

A. MME 

The algorithm first calculates the maximum and minimum 
eigenvalue of sample covariance matrix, namely, Xmax and 
\„^in■ If there is no signal, \nax/Kiin = 1; otherwise, 
^max/Knin > 1. So the resulting detector computes the ratio 
and is compared with a threshold 



T, 



MME 



> iMME 



(33) 



(23) B. AGM 



AGM is derived from GLRT, it finds an unstructured esti- 
mate of Rj; to be Rg + (7,^1 under "Hi and ct^I under T-Lq. 
Am represent all eigenvalues of sample covariance matrix and 
M is the dimension of sample covariance matrix. Since the 
arithmetic mean is larger than geometric mean, the resulting 
detector computes the AGM of the eigenvalues of sample 
covariance matrix and compares with a threshold 

TaGM — ''^^"^ — T- > lAGM (34) 



(n™A™)^ 



C. FTM 



FTM extracts leading eigenvector as the feature rji, which 
is stable for signals and random for noise. First, it learns the 
feature blindly by comparing similarity of two consecutive 
sensing segment, namely feature learning algorithm (FLA). 



FLA 



max 
;=i,2,---L-fe+i 



fe=i 



m [k] Vi+i [k + 1] 



(35) 



If Tpi^A > 7e feature is learned as (pg i = ?7,;+i. Then with the 
learned signal feature (ps^i as prior knowledge, this algorithm 
just compares the feature i from new sensing segment and 
signal feature to determine if signal exists. 



TpTM 



max 

1=1,2, ■■■L-k+l 



L 

E 

k=l 



<PsA [k] 4>xA [k + I] 



> IFTM 

(36) 



D. EC 

EC ifTOl method assumes the signal follows zero mean 
Gaussian distribution with covariance matrix Rs, and noise 
follows zero mean Gaussian distribution with covariance ma- 
trix (T^^I, 

s~A/'(Q,R,),n-A/'(G,o-2l) (37) 
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TABLE I 

Gain between FMD and EC using different number of data 



Number of samples 


2.5e4 


5e4 


le5 


2e5 


SNR of FMD 


-26 dB 


-30 dB 


-31.5 dB 


-34 dB 


SNR of EC 


-22 dB 


-23.5 dB 


-25 dB 


-26.5 dB 


Gain 


4 dB 


6.5 dB 


6.5 dB 


7.5 dB 



Both Rs and ct^ are given priorly. The hypothesis "Hi is true 
if: 

Tec = Yl ^J^- + ^"I) ' > ^I'c (38) 

E. Computational Complexity 

The computation complexity comes from two parts: sample 
covariance matrix calculation and eigenvalue decomposition of 
covariance matrix. Besides LNg multiplications and additions 
for MME, AGM, FTM, EC to calculate sample covariance 
matrix, FMD requires K additional additions to calculate 
covariance matrix. Generally O (L"^) multiplications and ad- 
ditions ai-e sufficient for MME, AGM, FTM to calculate 
eigenvalue decomposition. While FMD needs 1 function of 
matrix calculation, computational complexity depends on what 
kind of function is chosen. 

V. Simulation Results 

In the following, we will give some simulation results using 
DTV signal (field measurements) captured in Washington 
D.C. lUII- The data rate of the vestigial sideband (VSB) DTV 
signal is 10.762 MSamples/sec. The recorded DTV signals 
were sampled at receiver at 21.524476 MSamples/sec and 
down converted to a low central IF frequency of 5.381119 
MHz. The number of total samples contained in the segment 
is Ns = 100, 000 (corresponding to about 5ms sampling time), 
and they are divided into 166 sub-segments with 600 samples 
in each sub-segment. The smoothing factor L is chosen to 
be 32. Probability of false alarm is fixed with Pfa = 1%. 
The SNR of the received signal are unknown. In order to use 
the signals for simulating the algorithms at very low SNR, 
we need to add white Gaussian noise to obtain various SNR 
levels. 2,000 simulations are performed on each SNR level. 

The probabilities of detection varied by SNR for FMD 
compared with EC, AGM, FTM, MME are shown in Fig. [T] 
FMD method is better than all other methods introduced in 
Section [rV] which can achieve 0.5 probability of detection at 
SNR -31.5 dB. In addition. The performance can be improved 
further if we use more data samples. We can see from Table [l] 
the SNR of FMD can even achieve -34 dB at 0.5 probability 
of detection employing 2e5 data. Meanwhile, the tendency of 
gain between FMD and EC is apparent that it is increasing 
with more data samples. This is because more samples can 
make the numerical difference between two hypotheses more 
stable, hence the threshold can totally separate signals from 
two hypotheses. 
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Fig. 1. Probability of detection at 1% probability of false alarm 



VI. Conclusion 

In this paper we considered the spectrum sensing for single 
PU with single antenna, and FMD algorithm is proposed. 
Captured DTV signal is simulated with our proposed algorithm 
and other state-of-the-art algorithms, including MME, FTM, 
AGM. The simulation results showed our proposed algorithm 
can work lower than -30 dB SNR with limited data, which can 
even beat the performance of EC with perfect prior knowledge. 

The optimal number for each subsegment of 100,000 total 
data is 600. We need to find out how to optimally choose the 
number of data to calculate covariance matrix under different 
total sample data in the future. 
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